Development of a refinement criterion for adaptive mesh refinement in steam-assisted gravity drainage simulation

Steam-assisted gravity drainage (SAGD) is an enhanced oil recovery process for heavy oils and bitumens. Numerical simulations of this thermal process allow us to estimate the retrievable volume of oil and to evaluate the benefits of the project. As there exists a thin flow interface (compared to the reservoir dimensions), SAGD simulations are sensitive to the grid size. Thus, to obtain precise forecasts of oil production, very small-sized cells have to be used, which leads to prohibitive CPU times. To reduce these computation times, one can use an adaptive mesh refinement technique, which will only refine the grid in the interface area and use coarser cells outside. To this end, in this work, we introduce new refinement criteria, which are based on the work achieved in Kröner and Ohlberger (Math Comput 69(229):25–39, 2000) on a posteriori error estimators for finite volume schemes for hyperbolic equations. Through numerical experiments, we show that they enable us to decrease in a significant way the number of cells (and then CPU times) while maintaining a good accuracy in the results.     

Spatial variation of erosion in a small,glaciated basin in the Teton Range,Wyoming, based on detrital apatite (U‐Th)/He thermochronology

Evolution of mountain landscapes is controlled by dynamic interactions between erosional processes that vary in efficiency over altitudinal domains. Evaluation of spatial and temporal variations of individual erosion processes can augment our understanding of factors controlling relief and geomorphic development of alpine settings. This study tests the application of detrital apatite (U‐Th)/He thermochronology (AHe) to evaluate variable erosion in small, geologically complex catchments. Detrital grains from glacial and fluvial sediment in a single basin were dated and compared with a bedrock derived age‐elevation relationship to estimate spatial variation in erosion over different climate conditions in the Teton Range, Wyoming. Controls and pitfalls related to apatite quality and yield were fully evaluated to assess this technique. Probability density functions comparing detrital age distributions identify variations in erosional patterns between glacial and fluvial systems and provide insight into how glacial, fluvial, and hillslope processes interact. Similar age distributions representing erosion patterns during glacial and interglacial times suggest the basin may be approaching steady‐state. This also implies that glaciers are limited and no longer act as buzzsaws or produce relief. However, subtle differences in erosional efficiency do exist. The high frequency of apatite cooling ages from high altitudes represents either rapid denudation of peaks and ridges by mass wasting or an artifact of sample quality. A gap in detrital ages near the mean age, or mid‐altitude, indicates the fluvial system is presently transport limited by overwhelming talus deposits. This study confirms that sediment sources can be traced in small basins with detrital AHe dating. It also demonstrates that careful consideration of mineral yield and quality is required, and uniform erosion assumptions needed to extract basin thermal history from detrital ages are not always valid.     

Requirements for Topology in 3D GIS

Topology and its various benefits are well understood within the context of 2D Geographical Information Systems. However, requirements in three‐dimensional (3D) applications have yet to be defined, with factors such as lack of users’ familiarity with the potential of such systems impeding this process. In this paper, we identify and review a number of requirements for topology in 3D applications. The review utilises existing topological frameworks and data models as a starting point. Three key areas were studied for the purposes of requirements identification, namely existing 2D topological systems, requirements for visualisation in 3D and requirements for 3D analysis supported by topology. This was followed by analysis of application areas such as earth sciences and urban modelling which are traditionally associated with GIS, as well as others including medical, biological and chemical science. Requirements for topological functionality in 3D were then grouped and categorised. The paper concludes by suggesting that these requirements can be used as a basis for the implementation of topology in 3D. It is the aim of this review to serve as a focus for further discussion and identification of additional applications that would benefit from 3D topology.     

Kinetics and mechanism of natural fluorapatite dissolution at 25 °C and pH from 3 to 12

The dissolution rates of natural fluorapatite (FAP), Ca₁₀(PO₄)₆F₂, were measured at 25 °C in mixed-flow reactors as a function of pH from 3.0 to 11.7, and aqueous calcium, phosphorus, and fluoride concentration. After an initial preferential Ca and/or F release, stoichiometric Ca, P, and F release was observed. Measured FAP dissolution rates decrease with increasing pH at 3 ? pH ? 7, FAP dissolution rates are pH independent at 7 ? pH ? 10, and FAP dissolution rates again decrease with increasing pH at pH ? 10. Measured FAP dissolution rates are independent of aqueous Ca, P, and F concentration at pH ≈ 3 and pH ≈ 10.Apatite dissolution appears to be initiated by the relatively rapid removal from the near surface of F and the Ca located in the M1 sites, via proton for Ca exchange reactions. Dissolution rates are controlled by the destruction of this F and Ca depleted surface layer. The destruction of this layer is facilitated by the adsorption/penetration of protons into the surface at acidic conditions, and by surface hydration at neutral and basic conditions. Taking into account these two parallel mechanisms, measured fluorapatite forward dissolution rates can be accurately described using

     

Rock-magnetic properties of topsoils and urban dust from Morelia (>800,000 inhabitants), Mexico: Implications for anthropogenic pollution monitoring in Mexico’s medium size cities

In this work, we investigate the correlation between some magnetic parameters and the level of contamination by heavy metals in urban soils from Morelia city, western Mexico. The magnetic study was carried out on 98 urban soils samples belonging to distinct land uses. Most of analyzed samples contain ferrimagnetic minerals as the responsible for magnetization, most probably corresponding to the titanomagnetites/ titanomaghemites solid solutions. This is inferred from the susceptibility vs. temperature measurements and the isothermal remanent magnetization (IRM) experiments. These measurements also indicate that most of samples are almost completely saturated before 300 mT. Additionally, the S_-200 values (S_-200 = IRM_-200/ SIRM, where IRM_-200= Back-field of 200 mT after magnetic saturation) are between 0.7 and 1.0, characteristic of low coercivity magnetic minerals. The averaged saturation isothermal remanent magnetization (SIRM) curves can be used as an indicator of pollution level, as these curves show different saturation values according to the level of contamination by heavy metals: Cu, Ni, Cr and Sr. These associations of (titano)magnetite with heavy metals were observed by Scanning Electron Microscope revealing some complex aggregates rather than commonly detected spherules.     

Archeointensity investigation on pottery vestiges from Puertas de Rolón, Capacha culture: In search for affinity with other Mesoamerican pre-Hispanic cultures

We report the results of a rock magnetic and archeointensity study on pottery fragments recovered at the archaeological site of Puertas de Rolón, Capacha culture (Colima), that is considered one of the most important pre-Classic cultural developments of western Mesoamerica, Mexico. Standard rock magnetic experiments were carried out in order to investigate the reliability of the ancient field strength recovered from studied materials. Six out of the seven analyzed fragments gave reliable intensity determinations. The mean archeointensity value obtained in this study is 35.3 ± 3.4 μT. It is in good agreement with another contemporaneous Mesoamerican mean archeointensity of 34.3 ± 1.1 μT obtained for El Opeño. Statistically similar rock-magnetic parameters between ceramics from El Opeño and Puertas de Rolón potsherds appear to be related to clay-sources with similar geological settings and weathering processes. Alternatively, ancient trade or interexchange between western (Capacha Culture) and eastern Mesoamerica locations, as it has been proposed earlier, seem to be supported by results obtained in this study.     

Fractional crystallization of high-K arc magmas: biotite- versus amphibole-dominated fractionation series in the Dariv Igneous Complex,Western Mongolia

Many studies have documented hydrous fractionation of calc-alkaline basalts producing tonalitic, granodioritic, and granitic melts, but the origin of more alkaline arc sequences dominated by high-K monzonitic suites has not been thoroughly investigated. This study presents results from a combined field, petrologic, and whole-rock geochemical study of a paleo-arc alkaline fractionation sequence from the Dariv Range of the Mongolian Altaids. The Dariv Igneous Complex of Western Mongolia is composed of a complete, moderately hydrous, alkaline fractionation sequence ranging from phlogopite-bearing ultramafic and mafic cumulates to quartz–monzonites to late-stage felsic (63–75 wt% SiO₂) dikes. A volumetrically subordinate more hydrous, amphibole-dominated fractionation sequence is also present and comprises amphibole (±phlogopite) clinopyroxenites, gabbros, and diorites. We present 168 whole-rock analyses for the biotite- and amphibole-dominated series. First, we constrain the liquid line of descent (LLD) of a primitive, alkaline arc melt characterized by biotite as the dominant hydrous phase through a fractionation model that incorporates the stepwise subtraction of cumulates of a fixed composition. The modeled LLD reproduces the geochemical trends observed in the “liquid-like” intrusives of the biotite series (quartz–monzonites and felsic dikes) and follows the water-undersaturated albite–orthoclase cotectic (at 0.2–0.5 GPa). Second, as distinct biotite- and amphibole-dominated fractionation series are observed, we investigate the controls on high-temperature biotite versus amphibole crystallization from hydrous arc melts. Analysis of a compilation of hydrous experimental starting materials and high-Mg basalts saturated in biotite and/or amphibole suggests that the degree of K enrichment controls whether biotite will crystallize as an early high-T phase, whereas the degree of water saturation is the dominant control of amphibole crystallization. Therefore, if a melt has the appropriate major-element composition for early biotite and amphibole crystallization, as is true of the high-Mg basalts from the Dariv Igneous Complex, the relative proximity of these two phases to the liquidus depends on the H₂O concentration in the melt. Third, we compare the modeled high-K LLD and whole-rock geochemistry of the Dariv Igneous Complex to the more common calc-alkaline trend. Biotite and K-feldspar fractionation in the alkaline arc series results in the moderation of K₂O/Na₂O values and LILE concentrations with increasing SiO₂ as compared to the more common calc-alkaline series characterized by amphibole and plagioclase crystallization and strong increases in K₂O/Na₂O values. Lastly, we suggest that common calc-alkaline parental melts involve addition of a moderate pressure, sodic, fluid-dominated slab component while more alkaline primitive melts characterized by early biotite saturation involve the addition of a high-pressure potassic sediment melt.     

Phlogopite- and clinopyroxene-dominated fractional crystallization of an alkaline primitive melt: petrology and mineral chemistry of the Dariv Igneous Complex,Western Mongolia

We present field relationships, petrography, and mineral major and trace element data for the Neoproterozoic Dariv Igneous Complex of the Altaids of Western Mongolia. This unique complex of high-K plutonic rocks is composed of well-exposed, km-scale igneous intrusions of wehrlites, phlogopite wehrlites, apatite-bearing phlogopite clinopyroxenites, monzogabbros, monzodiorites, and clinopyroxene-bearing monzonites, all of which are intruded by late stage lamprophyric and aplitic dikes. The biotite-dominated igneous complex intrudes depleted harzburgitic serpentinite. The observed lithological variability and petrographic observations suggest that the plutonic rocks can be ascribed to a fractionation sequence defined by olivine + clinopyroxene ± Fe–Ti oxides → phlogopite + apatite → K-feldspar + plagioclase → amphibole + quartz. Notably, phlogopite is the dominant hydrous mafic mineral. Petrogenesis of the observed lithologies through a common fractionation sequence is supported by a gradual decrease in the Mg# [molar Mg/(Fe_total + Mg) × 100] of mafic minerals. Crystallization conditions are derived from experimental phase petrology and mineral chemistry. The most primitive ultramafic cumulates crystallized at ≤0.5 GPa and 1,210–1,100 °C and oxygen fugacity (fO₂) of +2–3 ?FMQ (log units above the fayalite–quartz–magnetite buffer). Trace element modeling using clinopyroxene and apatite rare earth element compositions indicates that the dominant mechanism of differentiation was fractional crystallization. The trace element composition of a parental melt was calculated from primitive clinopyroxene compositions and compares favorably with the compositions of syn-magmatic lamprophyres that crosscut the fractionation sequence. The parental melt composition is highly enriched in Th, U, large ion lithophile elements, and light rare earth elements and has a pronounced negative Nb–Ta depletion, suggestive of an alkaline primitive melt originating from a subduction-imprinted mantle. Comparison with a global compilation of primitive arc melts demonstrates that Dariv primitive melts are similar in composition to high-K primitive melts found in some continental arcs. Thus, the high-K fractionation sequence exposed in the Dariv Igneous Complex may be a previously unrecognized important fractionation sequence resulting in alkali-rich upper crustal granitoids in continental arc settings.     

Silicon in iron meteorite metal

Abstract– We report Si concentrations in the metal phases of iron meteorites. Analyses were performed by secondary ion mass spectrometry using a CAMECA 1270 ion probe. The Si concentrations are low (0.09–0.46 μg g^?1), with no apparent difference in concentration between magmatic and nonmagmatic iron meteorites. Coexisting kamacite and Ni‐rich metal phases have similar Si contents. Thermodynamic calculations show that Fe,Ni‐metal in equilibrium with silicate melts at temperatures where metal crystallizes should contain approximately 100 times more Si than found in iron meteorites in this work. The missing Si may either occur as tiny silicate inclusions in metal or it may have diffused as Si‐metal into surrounding silicates at low temperatures. In both cases, extensive low‐temperature diffusion of Si in metal is required. It is therefore concluded that low Si in iron meteorites is a result of subsolidus reactions during slow cooling.     

Ensemble properties of comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u, g, r, i, z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(<H) ∝ 10^{(0.49±0.05)H} for H < 18, with evidence of a much shallower fit N(<H) ∝ 10^{(0.19±0.03)H} for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g ? r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.